THE SPHERICAL COLLAPSE MODEL IN A UNIVERSE 
WITH COSMOLOGICAL CONSTANT 



o 

O ■ E. L. LOKAS 

' Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland 

E-mail: lokas@camk.edu.pl 

> 

Q ■ Y. HOFFMAN 

, Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel 
' ■ E-mail: hoffman@vms.huji.ac.il 

We generalize the spherical collapse model for the formation of dark matter halos 
to apply in a universe with arbitrary positive cosmological constant. We calculate 
the critical condition for collapse of an overdense region and give exact values of 
^ ■ the characteristic densities and redshifts of its evolution. 

in 

[ 1 Introduction 

The spherical cohapse model was first developed by Gunn & Gottlll for a flat 
' universe with no cosmological constant. It assumes that the process of forma- 

. tion of bound objects in the universe can be at first approximation described 

by evolution of an uniformly overdense spherical region in otherwise smooth 
^ ^. background. Despite its simplicity, the model has been widely used to ex- 

plain properties of a single bound object via extensions such as the spherical 
infall model as well as statistical properties of different classes of objects via 
^ ' Press-Schechter-like formalisms. Given the presently mounting evidence for 

Ci ■ cosmological constant we generalize the model here to include its effect. 
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•'^ ■ 2 The cosmological model 

XI ■ 

' The evolution of the scale factor a = R/Ro = 1/(1 + 2;) is governed by equation 



da Ho 



l + r!o --1 +Ao(a2-l) 



-1/2 



(1) 



The parameters Hq and Aq are the present values of fl and A = A/(3i/^) 
where A = const is the standard cosmological constant. The evolution of 
these parameters with redshift z is given by n{z) = flo{l + z)^[Ho/H{z)]'^ and 
A(z) = Xo[Ho/H{z)]^ where [H{z)/Ho]^ = r!o(l+z)3-(r!o+Ao-l)(l+z)2+Ao. 
The linear growth of density fluctuations is given by 

D{a) = -^ rf{a)da (2) 
2a/ (a) Jo 
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where /(a) was defined in Eq. (^. The expression in (g) is normahzcd so that 
for ^ — 1 and A = we have D{a) = a. For arbitrary parameters (r2oi-^o)j 
expansion of the right-hand side of Eq. (^) around a = gives 

D{a) = a + Oia^). (3) 



3 Evolution of the overdense region 

We assume that at some time corresponding to redshift Zj the region of proper 
radius rj is overdense by Aj — const with respect to the background, that is 
it encfoses a mass M(ri) — 47rpb,i?'f (1 + Ai)/3, where /9b,i is the background 
density at ti. 

Evolution of this region is governed by the famihar energy equation 
1 fdrV GM Ar2 



which, with a new variable s — r/ri, can be rewritten in the form 

i + ni (- -i) +Xi{s^ -1) 



ds Hi r /I \ 1-1/2 

— - = where g{s) 

dt g{s) 



(5) 



and we used the notation Hi = H{zi), fij = ^l{zi) and Ai = A(zi). 

Assuming conservation of energy we find that the maximum expansion 
radius rta (or equivalently, sta = ?'ta/n) of the overdense region must obey the 
condition 

6lS?a + 52Sta + &3 =0 (6) 

where bi = Ai, 62 = 1 ^ + Ai) — Ai and 63 = f2i(l + Ai). The only solution 
to Eq. which is real, positive and reproduces the Aq = case in the limit 
of small Ao is 



1/2 



27r 



with <j) = arccos[a;/(x2 + ?/2)i/2]^ ^ ^ -Qbl^'^bs and y = [3(-46i - 276i6i)]i/2, 
while for Aq = we simply get Sta = —b^li>2- 

The condition for the general solution (|^) to be real is 

Ai > Ai,er = j^P(^i) - 1 (8) 

where 

5Ai 3Ai(8 + Ai) ^ 3g(Ai) 



and 

(?(Ai) = {Ai[8 - A? + 20Ai + 8(1 - Xif/^]}'/'. (10) 

Ai^cr is the critical density for the overdense region to turn around. In the 
hmit of Aq ^ we have p{Xi) — > 1 and we reproduce the well known condition 

Ai > Ai,„(Ao = 0) = - 1- (11) 

After turn-around the proper size of the region r{t) evolves almost indepen- 
dently of the presence of the cosmological constant and can be well approx- 
imated by the analytical solutions for Aq — 0: r/rm = (1 — cos0)/2 and 
t/tcoii = {9- sin6l)/(27r) with Q<9 <2ti. 



4 The characteristic densities 

Integrating equations (|l|) and (||) we get 



/ f{a)Aa = Hot, f g{s)ds = Hit. 

Jo JQ 



(12) 



Eliminating t we obtain equations which can be used to calculate the scale 
factor at turn-around (ata) and collapse (ocou) of a region with particular Aj 
at a given Zi 



/(a)da = ^ I gis)ds (13) 

-"i Jo 

f[a)Aa^2^J\{s)ds (14) 

where Sta is given by Eq. (0). 

Assuming that the mass inside the overdense region does not change, the 
overdensity inside the sphere of size r with respect to the background density 
at any time is 

5=^-i-4f-)'(i+^i)-i- (15) 

At early times, t — s- 0, we can expand the expressions on the left-hand sides 
of both equations in ( p^ ) around a = and s = respectively. Integrating 
term by term we can express time as a power series of a and s respectively. 
Inverting those series we obtain a and s as power series of t^/^ . Inserting them 
into Eq. ( [l5| ) and replacing i-dependence with a we find 

J = /i(r!o,Ao,Ai,Zi)a-}-0(a2), (16) 
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Figure 1: Left panel: parameter <5c as a function of Qq for ^coll = 0- Solid lines correspond 
from bottom to top to Aq = 0, 0.2, 0.4, 0.6, 0.8 and 1. The dashed line shows results for the 
flat case Qo + Aq = 1. Right panel: parameter 5c as a function of ^coii for four models 
with flo = 0.3 and Ao = 0,0.4,0.7 and 1. The thin straight line marks the fiducial value 
(5c = 1.68647 for Oq = 1, Aq = 0. 



where 

/i(rio, Ao, Ai,zi) 



1 - r!o - Ao ^ [Qiil + Ai) + Ai - 1](1 + Zi) 



fin 



f7i(l + Ai)2/3 



Given the behavior of the hnear growth factor D{a) in Eq. 
density contrast as predicted by hnear theory 

(5l = /i(17o,Ao, Ai,Zi)D(a). 



(17) 



we obtain the 



(18) 



A particularly useful quantity is the linear density contrast at the moment 
of collapse i.e. when s reaches zero 



Sc = h[no, Ao, Ai(acoii), Zi]D{acoii)- 



(19) 



Ai(acou) in the above equation means that Ai corresponding to Ocoii has to be 
determined numerically for a given zi from the Eq. (14). Numerical results for 
i5c are shown in Figure ^ In the special case of fig + Aq = 1 we reproduce the 
results of Eke, Cole & FrenkH, while for open universes with no cosmological 
constant our results match those of Lacey & Coleu. Although it is not obvious 
from Figure |l|, in the limit of fio ^ we have i5c 3/2, independently of Ag. 
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Figure 2: Left panel: parameter Ac as a function of Qq for Zcoll = 0. Solid lines correspond 
from top to bottom to Aq = 0,0.2,0.4,0.6,0.8 and 1. The dashed line shows results for the 
flat case Oq + Aq = 1. Right panel: parameter Ac as a function of Zcoii for four models 
with Ho = 0.3 and Aq = 0, 0.4, 0.7 and 1. The thin straight line marks the fiducial value 
Ac = 177.653 for Oq = 1, Aq = 0. 



Another useful quantity is the ratio of the density in the object to the 
critical density at virialization 

Ac = (acoii) = — 3 [1 + Ai(acoii)J (20) 

Pcrit Scoll V / 

where Scoii = '"coii/n and TcoIi is the effective final radius of the collapsed object. 
We assume that the object virializes at icoU, the time corresponding to s — > 0. 
Application of the virial theorem in the presence of cosmological constant leads 
to the equation for the ratio of the final radius of the object to its turn-around 
radius F = rcoii/Ua obtained by Lahav et ala: 2r]F^ — (2 + r])F + 1=0 
where f] = 2AiSta/[^i(l + Ai)]. Numerical results for Ac in different models 
are shown in Figure |^. Again, we agree with the results for the special cases of 
f2o Ao = and fio + Aq = 1 derived previously by Lacey & Colecl and Eke 
et ala respectively. It should be emphasized that the value of Ac is strongly 
model-dependent, contrary to the common assumption of Ac ~ 200. 



5 The characteristic redshifts 

It is sometimes useful to be able to estimate the redshift of a particular stage 
of evolution of perturbation given its present overdensity as predicted by linear 
theory, 6h{a = 1) = So- For the redshift of collapse combining equations (18) 
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Figure 3: Redshift of turn-around (left panel) and collapse (right panel) of density fluctuation 
with present density contrast Sq for four models with Qg = 0.3 and Aq = 0,0.4,0.7 and 1. 
The thin straight line marks the fiducial case of f2o = 1, Ao = 0. 



and ( pj| ) gives Sq — Sc{acoii)D{a = l)/_D(acoii)- Using the previously obtained 
results for and formula we can calculate the present linear density 
contrast of fluctuation that collapsed at ZcoU- This relation can only be inverted 
analytically in the case of fig = 1, Aq = when we get Zcoii — Sq/Sc — 1. In 
an analogous way one can obtain the redshifts of turn- around, zta- Numerical 
results for both quantities are shown in Figure ^. Usefiil fitting formulae for 
those results will be presented in our forthcoming paperd. 
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